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I In a recent study we have obtained correction terms to the large N asymptotic expansions of 

Ps| ' the eigenvalue density for the Gaussian unitary and Laguerre unitary ensembles of random 

{^jQi N X N matrices, both in the bulk and at the soft edge of the spectrum. In the present 

^ , study these results are used to similarly analyze the eigenvalue density for Gaussian and 

Laguerre random matrix ensembles with orthogonal and symplectic symmetry. As in the 
case of unitary symmetry, a matching is exhibited between the asymptotic expansion of the 
bulk density, expanded about the edge, and the asymptotic expansion of the edge density, 
^-H . expanded into the bulk. In addition, aspects of the asymptotic expansion of the smoothed 

, density, which involves delta functions at the cndpoints of the support, arc interpreted mi- 

! croscopically. 
^ ; PACS numbers: 02.50.Cw,05.90.+m,02.30.Gp 

o ■ 

^ ; 1 Introduction 

Perhaps the best known result in random matrix theory is the Wigner semi-circle law (see 
"5 ■ 6.g. [3]). Consider a real symmetric matrix, with elements on the diagonal i.i.d. random variables 

g . having finite variance and similarly the elements above the diagonal. The Wigner semi-circle law 

• • . tells us that after appropriate scaling, the limiting eigenvalue density is given by the semi-circle 

■ functional form 

cq\ 1 0, |A| > 1. ^ ' 

As a concrete example, the Gaussian orthogonal ensemble (GOE) of real symmetric matrices is 
specified by its diagonal entries being distributed according to the normal distribution N[0, 1] 
and its upper triangular entries according to N[0, l/-v/2]. Let p^^\X) denote the eigenvalue 
density for x iV matrices from the GOE. After the scaling V2N p^^\\/2N X) ^ Np{X) the 
N ^ oo limiting form of p^^^X) is given by ()1.1() . 

The functional form (|1.1|1 has implications with respect to averaging a so called linear statistic 
A = X^jLio(Aj) over the spectrum of random real symmetric matrices. Thus, if the N ^ oo 
scaling is such that ^p^^\aNX) p^J^{X), a{aNX) — > a(A) for some on then 

/oo rl 
p^^\X)a{X) dX ^ N pw(A)a(A)dA. (1.2) 
-oo J-1 

The result H1.2() immediately draws attention to corrections to the Wigner semi-circle law. Indeed 
in studies of the full distribution of linear statistics averaged over the GOE and certain of its 
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generalizations, it is necessary to compute the 0(1) term in the asymptotic expansion of ()1.2() 
jni- For this one seeks the asymptotic expansion of p^^\X), where p^^\X) is the signed measure 
(smoothed density) such that 

p^^\X)a{X)dX= / p^^\X)a{X)dX 



to all orders in the corresponding asymptotic expansions. In the case of the GOE itself the 0(1) 
term is known |171 [21 IH] , and one has 



^p-W(V^A) ~ /'w(A) + ^Q(5(A-l) + 5(A + l))-i^^=l=X|A|<i) (1.3) 



where xt = 1 if T is true and xt = otherwise. 

The correction term in (|1.3() exhibits a most remarkable feature, namely delta functions at 
the edge of the support of the spectrum. The appearance of the delta functions at a microscopic 
level, when one seeks directly the asymptotic expansion of ^p(^)(^/2iVA) rather than the 
asymptotics of the smoothed quantity p^^-*(A) has not, to the best of our knowledge, been 
previously studied. One of the purposes of this paper is to undertake such a study for the 
Gaussian and Laguerre ensembles in random matrix theory. Each of the three symmetry classes, 
orthogonal (/3 = 1), unitary (/3 = 2) and symplectic (/3 = 4) will be considered. For the Gaussian 
ensemble it is known ^ that (|1.3)) then generalizes to read 

^.•"'(V/SVA) ~ pw(A) + l(i-i)(i(i(A-l) + *(A + l))-i^,„.,<,)^ (1.4) 

Our task is to relate this expansion to the asymptotic expansion of the density itself. 

The expansion ()1.4|) clearly shows both a bulk effect and an edge effect. This is in keeping 
with there being both a (global) bulk regime, and an edge regime which must be treated sepa- 
rately in the asymptotic analysis. As these expansions relate to the same quantity, one would 
expect there to be a matching in an appropriate limit. This topic, initiated in j8j for the GUE 
and LUE, is another main theme of the present work. 

We begin in Section 2 by recalling the results from ^ relating to the asymptotic expansions 
of the global bulk density, and the soft edge density, in the GUE and LUE. We then proceed 
to write down higher order terms in these asymptotic expansions (obtained from the method of 
(Hj). These higher order terms are then used to further investigate the matching phenomenon 
alluded to above. 

In Section 3 formulas required in the study of the asymptotics of the density in the Gaussian 
and Laguerre ensembles with orthogonal and symplectic symmetry are gathered. These formulas 
are used in Section 4 to study the corresponding global density asymptotic expansions, and in 
Section 5 to study the soft edge density asymptotic expansions. In Section 6 we use the results 
of Section 2, 4 and 5 to study our main topics of interest, namely the matching between the 
bulk and edge asymptotic expansions, and the microscopic origin of the delta functions in (|1.4() 
and its Laguerre analogue. 
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2 The Gaussian and Laguerre ensembles with unitary symmetry 



2.1 Definitions and summary of known results 

The Gaussian unitary ensemble consists of the set of Hermitian matrices with diagonal entries 
distributed according to the normal distribution N[0, l/\/2] and with upper triangular entries 
distributed according to N[0, 1/2] + iN[0, 1/2]. The corresponding eigenvalue p.d.f. is given by 

1 ^ 

cU^'^' n i^k-x,f (2.1) 

1=1 l<j<k<N 

where here and below C denotes some normalization constant. 

The Laguerre unitary ensemble can be specified by matrices X = A^A where A is a n x m 
(n > m) complex Gaussian matrix with entries distributed according to N[0, 1/ \/2] +^N[0, 1 / \/2] . 
All the eigenvalues are non-negative and have the joint distribution 

^ m 

-n^re-^' n i^k-x,)' (2.2) 

1=1 l<j<k<m 

where a = n — m. 

The eigenvalue p.d.f.s (|2.H) and (|2.2|) are special cases of the functional form 

1 ^ 

UE^(52) := ^n^^2(xO n (^k-xjf (2.3) 

1=1 l<j<k<N 

defining a general matrix ensemble with unitary symmetry in terms of its eigenvalue p.d.f. Thus 

Ev(GUEiv) = UEjv(e"^''), Ev(LUEjv) = UEiv(x°e-^') 

where Ev(M) denotes the eigenvalue p.d.f. of the matrix ensemble M. 

It is a basic result in random matrix theory (see e.g. fS!) that the eigenvalue density for 
the ensemble (j2.3p can be written in terms of the monic polynomials {pnix)} orthogonal with 
respect to the weight g2{x). Thus with {pn,Pn)2 denoting the corresponding normalizations we 
have 

p{x; UE;v(92)) = g2ix) (2.4) 

In a recent study Kalisch and Braak 10 have obtained the leading correction term to the 
Wigner semi-circle law for the asymptotic expansion of H2.4() in the case of the GUE. 

Proposition 1. Let — 1 < X < 1 be fixed. One has 

^ n(Y-JW (.-^N-'W n (X\ 2cOs(2iV7rPw(^)) 1 I \ .^ 

-p(X,UE^(e ))~Pw(X) ^3^^-^^^ N+^W^)^ (^-^^ 

where pw{x) given by li.ij) and 



X 1 

Pw{x) = 1 + —pwix) Arccosx. (2.6) 

2 vr 
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The methods in ^U] are particular to the Gaussian ensembles, relying on an integral formula 
coming from the supersymmetry method. Subsequently the present authors jSI have introduced 
a different strategy which reclaims ()2.5() . and furthermore applies equally well to the Laguerre 
case, for which the following result was obtained. 

Proposition 2. Let < X < 1 be fixed. We have 
lp(X;UE^(a;"e-^^^)) ~ pmp{X) 



cos{{2N + a)TTPMp{X)-aTT{l + XpMp{X)) a \ 1 / 1 v 

7r3X2p2^p(X) TT^XpMpiX)) \ } 



2.7) 



where 



2/1 2 

PMp{x) := -\ 1, Pmp{x) = I + xpMp{x) Arccos\/i (2. 

IT \ X vr 



(the subscript MP denotes Marcenko-Pastur, who first derived the limit law giving the leading 
term in this expression) . 

The strategy of |H] was to use integral representations of the product of orthogonal polyno- 
mials which result from the Christoffel-Darboux summation of (|2.4() . In addition to yielding the 
asymptotics of the global bulk density, it also gave asymptotics of the density at the so called 
soft edge. This is the name given to the boundary of the support at leading order, with the 
feature that the density, appropriately scaled, is non-zero on both sides of this boundary. 

Proposition 3. Let ^ be fixed. For the GUE 



20 

while for the LUE 



{^eikmf - 2i[k\'{i)f - 3Ai(e)Ai'(0) ]^ + O (^) ' (2-9) 



p 1+ . J,,,„ ;UEjv(3:"e- 



(2A^)2/3^V (2iV)2/3 
= ([Ai'(0]'-^[Ai(0]') + ^[Ai(e)]^ ^ 



Ari/3 

+^[2>e[m)f - ^m'm^ + (2 - ba^)M{iW{0) ]^ + ^ (^J • (2.10) 

2.2 Matching of the bulk and edge expansions 

We pursue the matching phenomenon observed in between the asymptotic expansion of the 
bulk density, expanded about the soft edge, and the asymptotic expansion of the edge density, 
expanded into the bulk. Explicitly, it was found that setting 

X = l+i/2N^I'\ (C<0) (2.11) 
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in the asymptotic expansion H2.5|) multiplied by N"^^^ /2, and then expanding in A^, a matching 
is obtained with the ^ — > — oo asymptotic expansion of the right hand side of (|2.9|) . This was 
checked on terms in ^ in the latter accessible by expanding the former; on this point we note that 
terms of all orders in inverse powers of in (\2.5\i . will after substitution of ()2.11|) contribute to 
each term in the expansion of (|2.9() . A similar matching was observed between the bulk LUE 
density expanded with the substitution (|2.11|) . and the ^ — > — oo asymptotic expansion of the 
0(1) and 0{N-^/^) terms of the edge density ((TTU)) . 

Based on this evidence, the hypothesis was put forward in [S] that the matching persists 
between all terms in ^, and at all orders in inverse (fractional) powers of N. Here we probe this 
hypothesis further by extending the asymptotic expansions (j2.5jl . (|2.7|) . (|2.9|) and (|2.1U|) . As 
explained in the orthogonal polynomial method readily allows for the computation of higher 
order terms, which we compute to be given as follows. 



Proposition 4. The 0{1/N'^) term in \2. 5\) is 



1 X{15 + 2X'^)sm[2N^TP^f^{X)]\ 1 

167r(l - A2)5/2 + 487r(l - A2)5/2 JJ^' 



the Oil/N^) term in lO) is 



the 0(1/A^2) term in fETj ) i 



IS 



a 



+- 



7r(l-X)2 

f8A 
192 

l + 4{-l + X)a'^ \ 1 



cos[2A^7rPMp(^) - 2QArccos VX] 



12 40 48/' ' " 40 ' ' ' ' V12 20 48/' ' " / A^ 
the 0{l/N) term in WUKl is 



a 



(2.12) 



180 + 981A2 + eOA^ + 4X6 ^ 

,3„,^(^_^,), c„sl2iV.Av(X)l^; (2.13) 

the oscillatory 0{\/N'^) term in \2. ,5)) is 

, 323190 + 647055X2 + 20358A'^ + 6084X6 - 1112A8 ^ 

-X ^ — — sin 2X7rPw(A) — r; (2.14) 

8294407r(l - X2 11/2 ^ wv ;j^4, v ; 



/-3 + 12X + 8X2 + 12(-l + X)(-l + 2X)a2N 
+ ( l92vr(l - Xfl^X^I^ ) ^-[2^-^MP(A) - 2aArccosyX] 



"647r(l-X)5/2x3/2jx2' ^^-^^^ 
the 0{l/N) term in is 

^ + ^ + fi) - ^Ai(OAi'(0 + f ^ - S - (Ai'(O)'l ^; (2.16) 



5 



Considering first the CUE, we now substitute (fTTT|) in (^3]) extended by ^TTI^ . (fTT^ . 
(|2.14|1 . Expanding the asymptotic series (an operation we denote by ~) gives the new asymptotic 
series in N, 



+ 16^|e|5/2 + ( " 2^ + 2304^0 cos(4|ep/V3) 48^1^^15/2^ ^) 

+ V 47r + 1287r|e|3/2 + I 4807r 331776^|epJ ^ ^ ^ 
233 |e|^sin(4|^|3/2/3)\ 1 



4608 20 / 7r|^|3/2 J 
O(^). (2.18) 



On the other hand, making use of the asymptotic series ^ 



— 2A:-1 

A;=0 k=0 



Aii-z) ~ 7r-i/2z-V4 ( sin(C + n/A) Y^{-lfc2kC'' - cos(C + n/A) Y^i'^)' C2k+iC 

(2.19) 



where 

S3, (216)'=/c! ^ - ^' 

the ^ ^ — 00 expansion of the 0(1) and ©(l/iV^/S) 

terms in (|2.9|) can readily be computed. 
Agreement is found with ah terms in (|2.18|) except the one involving the fraction 23695/331776. 
Thus, even though no terms 0{1/N) have yet appeared, whereas ()2.16() has a term at this 
order, the evidence is still in favor of a matching between all terms in ^, and at all orders in 
inverse powers of N. However this matching cannot be fully exhibited at any order in in the 
expansion of (|2.9|1 without knowledge of all terms in the asymptotic expansion of (|2.5|) . 

As already mentioned, a similar matching phenomenon was observed in |H] in the case of 
the LUE, and conjectured to hold at general orders as for the GUE. Further evidence for this 
conjecture can be obtained by substituting 1)2. 11(1 in 1)2. 7(1 extended by (|2.15|) . expanding as a 
series in inverse powers of N, and comparing against the ^ — > —00 expansion of (|2.1()|1 . The 
former operation gives 

' rp(l + ;rTL;UE^(."e-^^- 



jV2/3^V 2A^2/3' 

. / 2^/=^|g|V^ cos(4|g|3/2/3) N a(l + sin(4|^|3/2/3)) 
V ^ 24/3^1^1 ; + 22/37r|^|V2Ari/3 

+ 160vr|g|3/2Ar2/3 " 20a^ + 801^1^ + 40(-l + 2a'm'/' cosim'^'^ 

+(5-20a2 + 16|eP)sin(4|e|3/V3)) +0(1). (2.20) 

The latter operation gives agreement with this expansion at 0(1), 0(l/A^^/3) ^^d for the terms 
involving factors of \(,\'^ at 0(l/A^^/3^_ This is consistent with the matching hypothesis. 
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3 The Gaussian and Laguerre ensembles with orthogonal and 
symplectic symmetry — general formulas 



The Gaussian orthogonal ensemble has already been defined in the Introduction. At the level 
of an eigenvalue p.d.f., the GOE can be defined by the joint distribution 

1 ^ 

1=1 l<j<k<N 

Likewise the Gaussian symplectic ensemble, Laguerre orthogonal ensemble and Laguerre sym- 
plectic ensemble can be specified either in terms of the distribution of certain classes of random 
matrices (real symmetric matrices in the cases of orthogonal symmetry, and Hermitian matrices 
with real quaternion elements in the cases of symplectic symmetry) , or in terms of the functional 
form of the eigenvalue p.d.f. (see e.g. ,5. ). Here we will note only the latter, which in the case 
of the GSE reads 

1 ^ 



c 

1=1 l<j<k<N 



for the LOE reads 



and for the LSE reads 



1 ^ 



6 \xi^ "^ili 



1=1 ^<j<k<N 
N 



1=1 l<j<k<N 

For the Laguerre ensembles one requires the eigenvalues be positive and thus x/ > (/ = 
1 , . . . , N) . Thus we see that if we define a matrix ensemble with orthogonal and symplectic 
symmetry by the eigenvalue p.d.f.s 

N 



OENigi) := ^Yi^ii^i) Yi \xk-Xj\ 



C 

1=1 l<j<k<N 



and 

TV 



SENigi) := ^Y{gA{xi) Y[ (^k 



C 

1=1 l<j<k<N 

respectively, then we have 

Ev(G0E7v) = 0Ejv(e-^'/2)^ Ev(LOEjv) = OE^(x°/2e-^/2) 

Ev(GSE7v) = SE7v(e-2^'), Ev(LSE^) = SEjv(x2"e-^'). 

In the work P, the eigenvalue densities for the ensembles OE(gi) and SE((74) were computed 
for all the so called classical weights 

e~^^/^, Gaussian 
^(a-i)/2g-a;/2 (^x > 0), Laguerre 



gi{x) 



(l_^)(a-i)/2(i + a;){6-i)/2(_i<^<l)^ Jacobi ^^'^^ 
(1 +x2)-(°+i)/2, Cauchy 
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(3.3) 



Hermite 
Laguerre 
(1 + x)''+\ Jacobi 
+ Cauchy 

in terms of a formula depending on the symmetry class only. Thus for ensembles OEtvC^i) with 
classical weights (|3.2j) one has 



p{x; OEN{gi{x))) = p{x; \JEn-i{92{x))) 

+7 w ^5i(a;)pjv-i(2;)- 

[PN~2,PN-2)2(PN-l,PN-l)2 2, 



sgn(x - t)pAr„2 {t)gi (t) dt, (3.4) 



while for ensembles SEjv(fl'4) with classical weights (|3.3p one has 

1 



p{x;S'EN{gi{x))) 



In (1231) and 



-p{x;\]E2N{g2)) 



-gi{x) 



C2N-1P2n{x) 



g2{x) 



'^{P2N, P2N )2{P2N-1, P2N-1 ) 2 
2 

, Hermite 
x°e~^, Laguerre 
{l-xY{l + xf, Jacobi 
(l + x^)-", Cauchy 



gi{t)p2N-i{t)dt. (3.5) 



(3.6) 



while 



1 , Hermite 
i, Laguerre. 



(3.7) 



(Pi,Pi)2 [ ■ 

The quantities {p„(x)} are the monic classical orthogonal polynomials with respect to the weights 
(|3.6|) . and {pn,Pn)2 the corresponding normalizations. Thus in the Gaussian case 

Pn{x) = 2-"if„(x), {Pn,Pn)2 = l^^l'^T'^uX (3.8) 

while in the Laguerre case 

p„(x) = (-l)"n!L^(x), (p„,p„)2 =r(n + l)r(a + n + l). (3.9) 

Essential tools in our subsequent analysis of the asymptotic forms of 1)3. 4() and ()3.5|) are 
particular asymptotic formulas for the Hermite and Laguerre polynomials. Consider first the 
bulk region. In the case of the Hermite polynomials, the formula is due to Plancherel and Rotach 
[13] . It tells us that with 



a; = (2n + 1 



,1/2 



COS ( 



e < (/) < TT — e, 



we have 



e-^'/2^„(x) = 2"/2+i/4(^!)i/2(^^)-i/4^g-^^)-i/2j^g.^(|(^^/2+l/4)(sin2</.-2(^)+37r/4)+0(n-i; 



Setting 

V2NX = {2{N + m) + 1)^/^ cos (p 
with —1<X<1 fixed we deduce from this that for m fixed 

/ O \ 1/4 om/2+7V/2 , , , / 1 \ \ 

H^^UV2NX) = Q ^^3^iV"/^-^/^(iV!)^/^e^^^i^i(X)(l + 0(l)) (3.10) 
where 

gl^^ljix) := cos (^Nxy^l - x'^ + {N + l/2)Arcsinx - N7r/2 -mArccosxj. (3.11) 

The Plancherel-Rotach formula (|3.1()|) was extended by MoeckUn to the Laguerre polynomials 
HU. With 

X = (4n + 2a + 2) cos^ e < < 7r/2 - en"^/^ 



it reads 



e-^/2^(f)(^) ^ (-l)"(7rsin</))-V\-°/2-i/4„-/2-i/4 

sin ((n + (a + l)/2)(sin 20 - 2(/>) + 3ir/4) + (nx)~^/20(l) 



Setting 

4nX = (4(n + m) + 2a + 2) cos^ </) 
with m fixed and e/n < X < 1 we deduce from this that 

^./2^-./2^(a) = (-l)«+™(2vrVX(l - X))-V2„«/2-l/2/^W (X) + 0(i)) (3.12) 
x=4nX \ ' n / 

where 

ff(^)„(X) := sin (2n(v^X(l - X) - Arccos \/X) - (2m + a + l)Arccos ^fx + 37r/4) . (3.13) 



We turn our attention now to the (soft) spectrum edge. In the Hermite case, the formula of 
Plancherel and Rotach tells us that with 

X = (2iV)i/2 ^ 2-1/2^-1/6^ (3_14) 

we have 

exp{-x^/2)HN{x) = ^^/^2^/2+V4(Ar!)i/2Ar-i/i2|Ai(t) + 0(Ar-2/3)|^ 
where Ai(t) denotes the Airy function. It follows from this that with x again related to t by 

eM-xVmN+^m{x) = (2X)™/2^i/42^/2+i/4(Ar!)V2^-i/i2|Ai(t) _ ^Ai'{t) + 0{N-^/')]. 

^ (3.15) 

In the Laguerre case, Szego gives that with 

x = 4X + 2a + 2 + 2(2iV)^/^t (3.16) 



we have 

e-^/2_^^(x) = (-l)^2-"-i/3iV-i/3j^Ai(t) +0(Ar-2/3)y 
It then follows that for fixed p 

e-^/^L^ , (x) = (-i)^+P2-°-V3Ar-i/3 ( Ai(^) _ + " + ^ Ai^(e) + 0(iV-^/^)) . 

(3.17) 

4 Asymptotic expansions in the bulk 
4.1 GOE and GSE in the bulk 

As with the GUE, the asymptotic expansion of the bulk density for the GOE and GSE has been 
carried out by Kalisch and Braak JOl- But as their method is particular to the Gaussian case, 
we give an alternative method which can be extended to the Laguerre case. 

Consider first the GOE. Substituting the appropriate formula from (|3.7() . and H3.8|) . in (|3.4() 
we see after minor manipulation that 

p{x- 0E^(e--'/2)) = p{x- UE^._i(e--^)) + ^^1^^21^-1)1 [ ^''^"^^-^^1) dt. (4.1) 
Also, making use of 1)2. 4|) shows 

2-(JV-i)g-x2 



p(x;UE^_i(e-^ )) = p(x;UEjv(e-^ )) - ^^^^^^ _ ^^^ {Hm^i{x)Y . (4.2) 

From ()4.ip and 1)4. 2() the following asymptotic formula for the bulk eigenvalue density is obtained. 
Proposition 5. Let —1 < X <1 he fixed. We have 

Proof. First we note that 

p{X- OE;v(e-^"')) = V2Np{V2NX- OE;v(e-"'/')), (4.4) 

and we proceed to analyze the large form of the right hand side using (|4.1j) . In relation to 
the latter, by a simple change of variables. 



/ e-''^^HN^2{t)dt = V2N e-*'/2^7V-2(t) 

Jo Jo 

while making use of ()3.1Up shows 

Vm e-''/^HM-2it) 
Jo 



dT 

t=V2NT 
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The leading contribution to the integral for large comes from the neighborhood of the end- 
points T = and T = X. About T = 

i'^im ~ cos(2iVr-(iV + m)7r/2 + 0(r2)), 

while about T = X 

9m,Ni'^) ~ cos (^NXVl -X^ + {N + l/2)ArcsinX 

-Ntt/2 - mArccos X + 2Ny^l - X^{T - X) + 0{{T - Xf)y 

Thus we have 

/■^ g%iT) 1 AH) t,.. 

I (l_r2)i/4 ~ 2Ar(l-X2)3/4^-2,^W (4.6) 

where 

9mN(^) '■= sin (^A^xa/i - ^2 + (A^ + l/2)Arcsin x - A^7r/2 - mArccosx^ (4.7) 

and use has been made of the fact that A^ is assumed even in (|3.4p . 
We read off from (IXTU]) that 

^^"*^"'""%-/^(A'0'/^«i';Vm(l + o(i)). (4.8) 



Making use of this together with (|4.5I) . (|4.6|1 and Stirling's formula we deduce 



2iv~i^i/2(^_2)! 7g 



/ e'' '^HN-2{t)dt 

Jo 



1 l%ixrg%iX)^^^.i_ 



7tV2N (1-^2) V VA^ 
1 9%{X) 



{Xgi%iX) + VT^g%iX)) (l + 0(^)) (4.9) 



7r\/2iV (1-^2) 

where the equality follows from simple trigonometric identities. 
For the second term in 1)4. 2(1 . use of ()4.8p shows 



;iHN-iix)) 



x=V2NX 



7rV2(iV-i)! 

And for the first term in (|4.2|) we know from 1)2. 5|1 that 

piVmx-UENie-^')) 

TT ^ Va^ 27r(l-X2) +^lvAf3/2y'-^*-''^ 

But 

cos(2A^A:Vi - A:2 + 2Af Arcsln AT - A^^) = 71^^(2 (ffi% (X))^ - 1) - 2X5i^J^(X)5i^J^ (AT) . 

'(4.12) 
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Substituting (|I?T^ in then adding this, and jnUI), and recalhng (g31) gives (|0)) . 

□ 

The result (|4.3() agrees with that computed by Kahsch and Braak in and is also consistent 
with H1.3() . We remark that in ^U] the 0(l/iV^) term is also given, being equal to 

3 + cos((2iV - 1) Arcsin X + 2NXVT^X^ 



167r(l - X2)5/2Ar2 8^(1 - X2)5/2Ar2 

In principle the present method offers a systematic approach to all correction terms. For this 
we need the explicit form of the higher order terms in (|3.1U|) . and these can in fact be calculated 
from the results in However we have not pursued such calculations. We remark too that 
a calculation of the non-oscillatory portion of 1)4. 13(1 is undertaken in j^; however the result 
obtained does not agree with (|4.13|) . 

We turn our attention now to the GSE. First we note that 

p{X; SE^(e-^^"')) = 2VNp{2VnX; SEjv(e-"')). (4.14) 

Regarding the right hand side, making use of (|3.7|1 . 1)3. 8() as well as the integral evaluation (see 
e.g. P) 

(4.15) 

gives 

p(x;SE;v(e-^')) = ^p(x; UE2;v(e-"')) 

-4vrV2(2Ar_i)!^V2 22A^-i(Ar-l/2)!-' J, ^ ' H^N ^^^^t) dt y (4.16) 

The asymptotic form of (|4.16l) can be calculated according to the strategy of the proof of 
Proposition IS] to give the following result for the bulk scaled density in the GSE. 

Proposition 6. Let —1 < X < 1 be fixed, and 5q^]y(X) be given according to VJ.ll]) . We have 

1 ,^ , _4Arx^A 2 r — ( 1 (-1)^\ 5o,27v(^) 

-p(X;SEMe )) - -Vl-^^-(^ + ^) (,_^.)V4 

+4;Jiv7r^ + ''(^)- ^'-^'^ 

Proof. Analogous to (|4.5p . it follows from 1)3. 1U() that 
while proceeding as in the derivation of 1)4. 6|) shows 

^ S%Nit) 1 1 (H) (^.tDl 

g (l_t2)l/4^^ 4iV(l-X2)3/4^-l-2ivl^J 4^ • 
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Thus, after making use too of Stirling's formula, 



/2N V^V4A^(l -X2)3/4^-i.2Wv / 4^ 
Since H3.1U|) gives 

we thus have that with x = the final line in (|4.1()|) has the asymptotic behavior 



4V^(1-X2)l/4V' ^^V(l-X2)3/4 

Further, we see from (|2.5|) and that 

^/9(2^/iVX;UE27v(e-"')) 

iV.^/p-^ ^/r^(2(ggj,(X))2-l) + 2Xggi,(X)ffgj,(X) ^ 1 



vr 87r^/]V(l - X2) ViV3/2, 

Adding these last two results, and recalling (|4.16l) and (|4.14() gives the stated formula. □ 

In ^U], at 0{\/N) only the non-oscillatory term is reported. We note too that the non- 
oscillatory term at 0{1/N) in (|4.17j) is consistent with (|1.4jl in the case /3 = 4. 

4.2 The LOE and LSE in the bulk 

We now apply the same strategy to the Laguerre case. For the LOE, substituting the appropriate 
formula from p.7p . and 1)3. 9p . into (|3.4() shows 



L%_2{t)t(^~'^/^e-'/^dt-2j Lf,^,{t)t(^-'y'e-'/'dt) (4.18) 

while 

Kx;UE^_i(x"e--)) = /.(x; UE^(x-e--)) - (x^/^e-^/^L^^iCx))'. (4.19) 

Proposition 7. Let < X < 1. We have 

lp(X; OE^(x("-i)/2e-2^-)) ~ pmp(X) - .^V"^, + (4-20) 
A' 27rA y^X(l - A) 
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Proof. By a change of variables 

p{X; OEjv(x("-^)/2g-27Va;)) ^ 4Np{4NX; OEiv(x("-i)/2e-^/2))^ (4.2i) 

so the task is to analyze the — > oo asymptotics of the right hand side of (|4.18|) with x replaced 
by mX. 

We know from ^1 that 

^ ' 2"/2-3/2r(Ar)r((a + iV)/2) ^ ' 

Regarding the second integral in H4.18() . we first write 

L%__2{t)t(--'^/^e-'/^dt=(^J +j )L%_^{t)&-^^l^e-'/^dt 
where < e ^ 1. In relation to the region [4e,4A^X], (|l'{.12)) tells us that for even 

L%_,itrl'e~'l\^^^^ = (2vr7r(r3Y))-V2Ar(-i)/2(^W^(T) + 0(^))^ 
Substituting in the integral and changing variables T = shows 

Jie J^/7/N \fs{l- S^yi'^ 

(4.23) 

In relation to the interval t G [0, 4e], we know that for N ^ oo 
where Jni^) denotes the Bessel function, and thus 

/•4e fie Jj. 

/ L^_2(t)i(''-^)/2g-t/2^^ ^ / J^(2(iVt)l/2) ^. (4.24) 

Jo Jo Vt 

We expect the leading contributions to come from the neighborhood of the upper terminal 
s = \fX in ()4.23p , and the lower terminal t = in ()4.24p (the integrands should connect smoothly 
from the lower terminal of (|4.23j) to the upper terminal of (|4.24|) ). Since about s = \fX 

5-^2iv(^^) ~ sin (2iV(V^(l - X)-Arcos \/x)-(a-3)Arcos \/x+3^/4+47VVr^^(s- Vx) 



we have 

} V5(l-s2)l/4 4iV7r^XV4(l_X)l/4 

where 

9rn^N(^) ■~ ^2A^( Y^x(l — x) — Arcos "v/x) — (a + 1 + 2m)Arcos \/x + 37r/4^ . 

Also 



/•CO 

/ Jait)dt = l, 

Jo 
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so we have 

Reading off the asymptotic form of the factor x'^^'^e~^^'^L'^_-^^{x)\x=4NX in (|4.18|) from we 
deduce from this that 

p(4iVX;OE^(.(-i)/2e--/^)) ~ UE^-i(-"e--)) " ^^^^^^^^ 

It remains to determine the asymptotic form of (|4.19l) . For this we use the analogue of ()4.21() , 
and (|l'{.12|) . to obtain 



p(4iVX;UE;v-i(x"e--)) ~ UE^(x"e-^^-)) - x^- X) ^^-''""™' ' 

Noting from the definitions and by a simple trigonometric identity that 

9%{X) = {2X - irg%{X) - 2V(1 - X)Xg%{X) 
= {2X - irg%{X) + 2V(1 - X)Xg%iX) 

we therefore have 

p(4iVX;OE^(x("-)/2e-/2)) ~ i^P(^; UE^(-'^e-^^^)) " ^^^^^^^^^^ (4-25) 



Now, with 



we have that 



AN,a{X) := 2N{y/X{l- X) - Arccos Vx) - aArccos Vx (4.26) 



g^^l^iX) = sin(A^,«(X) + Arccos + 3^/4) 
g^^^l{X) = cos{AN,aiX) - Arccos Vx + Sir/A) 



and thus 



<7%(X)5S(^) = ( - cos(2A^,„(X)) + 2VX(1-X)) . (4.27) 



1 
2 

Substituting (fTTTI) in and noting from (|T7|) that 

— p(X; UE^(x"e-^^-)) ~ ipMP(X) - ""^^^^^'"^^^ + " (4.28) 

4Af^^ ' " 4^^^^ ^ 167rX(l - X)N ^ Sir^Xil - X)N ^ ' 

we obtain □ 

Consider now the LSE. Analogous to ()4.21|1 . by a change of variables 

y9(X;SE7v(x°+^e-*^^)) = 87V/5(87VX; SE7v(x°+^e-^)), 

while (pT3|) together with the fact US H] 

/>oo 

/ e-*/¥"-i)/2^f^_^(^)^^^0 (4.29) 

JO 
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shows that 



x/2^(a-i)/2^a^(^) / e-*/2t("-i)/2L^^_^(t)dt. (4.30) 



r(i + 2iv) 

~4r(a + 2iV)^ 

Proposition 8. Let < X < 1. In terms of the notation iS.lclp and \4-^^\ ) '^^ have 
^p(X;SEAr(:r"+ie-S^^)) ~ pMp(X) 

' + " + ^ +o(iV-^). (4.31) 



2(7riV)i/2x3/4(l_X)V4 47riVV^(l-X) 
Proof. Following the strategy of the proof of Proposition [7| we find 

io ' * L^^-iW^t ~ (2iV)^ + i2^y/2 2iVXV4(i_x)3/4- 

Also, (|3.12|) and Stirling's formula show 

1 &{X) 



r(i + 2iv)_^_,/2^(„_,)/2^„^(^) 



4r(a + 2iV) ^^''x=87vx 16(7riV)i/2(2Ar)(°-i)/2x3/4(l-X)V4- 

Substituting in (|4.3Up gives 

p{8NX-SEN{x^+'e--)) ~ -J_p(X; UE2jv(x"e-8^^)) 



16(^iV)i/2x3/4(l-X)i/4 32^iVX(l-X) 
Analogous to 1)4. 27() we have 

&iXr9%N{X) = -1(cos2A2nAX) + 2y/X{l-X)). 



2 

Substituting this together with in (|0^ gives (|OT|) . □ 



5 Asymptotic expansions at the edge 
5.1 The GOE and GSE 

The scaled densities p{X;OE]\i{e~'^'^^)) and p{X;SE]\[{e~^^-^^)) have to leading order their 
support on (—1, 1). We know from previous studies [HI El that setting X as specified by 1)2. 11(1 
(with the restriction ^ < removed) and multiplying by N^^^, the limit — > oo exists and can 
be computed explicitly. Here we are interested in computing the first correction, as in the soft 
edge formula (|^ for the GUE. 

In the case of the GOE, we see from (|4.4j) . (|4.H) and (|4.2j) that in addition to the knowledge 
of (E31), an asymptotic formula for p(l + ^/2A^2/3. 0E7v(e-^''')) can be obtained by making use 

of (Trr^ . 
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Proposition 9. We have 



2Af2/3 



Proof. Consider the integral in ()4.1() . We know from 1^ that 

Replacing by A^ — 2, setting x = (2A^)^/2 + 2^^/2A^^-'^/^^, making use of ()3.15() and simplifying 
shows that 

»(2Ar)i/2+2-i/2Ar-i/6j 

Jo 

Now using p.lSjl with m = — 1, and multiplying with the result H5.2|l as required by (|4.1j) we 
obtain 



2^-%V2(Ar_2)! 7o 



'HN-2{t)dt 



2;=(2Ar)i/2+2-i/2Ar-i/6g 



"^'^^ 'Ai(0 + ^Ai'(o) (l - r Ai(2/) dy + ^Ai(e)) . 



2V2 

According to H4.2|) . we also require the asymptotic formula 

(e-^V2^^_,(x))^ ~ -^(Ai'(0)^ 



7rV2(Ar_l)!^ a;=(2Af)i/2+2-i/2Ar-l/6g JS[l/& 

which follows from H3.15() . Further 

(cf. ()4.4() ) so the asymptotic form of the first term on the the right hand side of 1)4.21) can be 
read off from (|2.9|) . Doing this, and recalling ()4.4p . we deduce (|5.ip . □ 

We turn our attention now to the GSE. We will analyze (|4.16p rewritten to read 
Proposition 10. We have 

~ (Ai '(0)' - C(Ai(0)' - ^Ai(e) ( Ai(t) dt - j^^m) + 0(Ar-2/3)^ . 

(5.4) 
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Proof. By a simple change of variables 



p 1 + 



2(2iV)2/3 



(5.5) 



so we seek the large asymptotic form of (|XT3|) with x ^ (4iV)^/2 _^ 2~^/'^{2Ny^l^^. Making 
use of ()3.15p shows 



Arl/6 , /-oo 1 , 



r={4Af)l/2+2-l/2(2Ar)-l/6^ 



(5.6) 



while it follows from (|2.9j) that 



lp((4iV)V2 + 2-V2(2Ar)~i/6^. UE2;v(e-^^)) ~ ((Ai'(O)^ - ^(Ai(0)^ + 0^-^^ 

Substituting ()5.7() and (|5.6() in 1)5. 3p and recalling 1)5. 5j) gives the stated result. 



(5.7) 

□ 



5.2 The LOE and LSE 

The soft edge scaling variables for the LOE and LSE are the same as those for the LUE, exhibited 
in 1)2.10(1 . The leading correction term to the corresponding soft edge densities are given by the 
following result. 



Proposition 11. We have 



-;OE^(x("-i)/2e-2^^; 



(2A^)2/3'^V (2iV)2/3 

~ (Ai'(O)' - miO? + ^Ai(0(l - Ai(s) ds 
{a -I) 



2(2A^)i/3 



Ai'(0(l- / A\{s)ds]-{Kmf +0(iV-2/3) 



(5.8) 



and 



9 , £■ , -1 poo 

p(l + 77^;SEjv(x-+ie-«^-)) ~ (Ai'(0)'-e(Ai(e))'--Ai(0 j Ai{s)ds 



(4Ar)2/3'^V (4iV)2/3 

+^^^ {mo? + Ai'(0 Ai(t) dt) + 0{N-^/'). 

Proof. Consider first the LOE. We rewrite (|4.18() as 



(5.9) 



p{x;OEN{x^--'y\^^/')) = p(x;UE^_i(x"e-^) + ^ 



, poo roo 

x(2j L%_^{t)t^"-'y'e-'/'dt- J L^_2(t)t("-i)/2e-*/2dt 



(5.10) 
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The asymptotic form of the final integral in H5.9() is known from 1)4. 22() . According to H4.21() , we 
seek the asymptotic form of the remaining terms with 

For the first integral, use of ()3.17|) gives 

4iV+2(2Af)l/3g 



iV(-i)/2 ( P Ai(.) ds - (|^Ai(0 + 0(iV-2/3; 



The asymptotic form of the term outside the bracketed integrals is given by (|3T7|l withp = -1. 
For the first term on the right hand side of H5.1U() . use of ()4.19|) . the analogue of ()4.21|1 . and 

shows 

The asymptotic form of all terms have now been determined, and 1)5.8(1 follows. 
In relation to the LSE, we use (|4.29p to rewrite (|4.3U() to read 



oo 



+^^^^"'^'"^""'^^'^^^(") / e-"'t^''-'^"L^^-S)dt- (5.11) 



A simple change of variables gives 



p(l + SE;v(x"+ie-«^-')) = mp[m + 2{ANf'H; SE^(x-+ie--') ) , (5.12) 

so we seek the asymptotic form of (|5.1ip with 

X = m + 2{mf/^£^. 
For this, we make use of (|3.17() to deduce 



4r(a + 2N) 

Further, making use of the analogue of (|5.12|1 and recalling ((2.1U|) shows 
^(SN + 2(4Af)'"{;UE2N(:c"e-)) 

The asymptotic form of all terms in ((5. 11(1 are now determined. After use of (|5.12() . (|5.9() follows. 

□ 

Figure 1 provides a numerical comparison of the asymptotic expansion given by ()5.9() with 
the exact result for the LSE density given by ()5.1H) and (j5.12|) . 
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Figure 1: Numerical comparison of the asymptotic expansion (|5.9|) . shown as the dashed hne, 
and the exact result for the LSE density given by H5.11() and H5.12|) . shown as the solid line, for 
the eigenvalue density near the soft edge of the LSE with = 20 and a = 1/2. 



6 Consequences 

6.1 Matching the bulk and edge expansions 
The GOE and GSE 

In Section r2.2l the conjectured matching between the bulk asymptotic expansion expanded about 
the soft edge, and the asymptotic expansion of the edge density expanded into the bulk, for 
both the GUE and LUE was given further credence by its confirmation on further terms in 
the asymptotic expansion. To study the matching in the case of orthogonal and symplectic 
symmetry, it is therefore convenient to express the corresponding asymptotic expansions in 
terms of the unitary symmetry expansions. We will consider first the Gaussian cases. 
According to (gSJ, (1031) . (|2T2|l 



-p{X;OEM{e 



+ 



)) ~ lp(X;UE^(e-^-^)) + 

1 + 2X2 x{2l + 2X2) sin[2iV7rPw(^)] 
87r(l -X2)5/2 487r(l - X^f/^ 



1 



cos[2iV7rPw(^)] 



1 

iV 



cos[2iV7rPw(^)] \ 1 , . 
87r(l -X2)2 y'iV2' ^ ' 



20 



according to ()5.1() and (|2.9|1 



2iV2/3'^V 2iV2/3' V 2iV2/3^V 2iV2/3' 



Ai(t)<it)+^Ai'(0(l-y^ Ai(t)dt)+0(^); (6.2) 
according to (|4.17|) and (|2.5jl 

SE;v(e )) ~ — VEMe )) + 3,(1 _ ^2)^ 
1 ^ 1 ^ 1 x cos[2iV^Pw(X) + iArcsinX] , ^(..-3/2^. (g 3^ 



AttNVI - X2 V^2^ 27riV; (1 - X2)i/4 

according to ()5.4() and (|2.9|) 



(2iv)2/3^V 2(2Af)2/3' V 2(2iV)2/3'^V 2(2iV)2/3 

-iAi(e) rAiW^t + T^^S + ofTTL); (6-4) 



2(2A^)V3 ViV2/3 
Substituting (|2.11|1 in (|6.1|) and and expanding as in (|2.18|) gives 

' rp(l + -L;OE^(e-^^^)) ~ J^pfl + :T4?;UE^(e-2^^-^: 



27V2/3^V 2Af2/3' V 2iV2/3^V 2iV2/3 



and 



3 23sin(4|ep/V3) cos(4|ep/V3) 

^167r|<ep/2 ^ 967r|ep/2 ^ 87r|^| 
/ 1 cos(4|^|3/2/3) N 1 

v4^ie|i/2 i6^ie|2 /Ari/3' ^"■''^ 



P{ 1 + TTTTTl^; SEjv(e-^^^" ) ~ TT^TTTT^P 1 + UE2^(e 



(2iV)2/3^V 2(2A^)2/3' V 2(2iV)2/3^V 2(2iV)2/3 

1 1 sin(||g|3/2+vr/4) cos(||ep/^) 
2^ |e|V4 + 87t\C\ 

+ f ^ + ^_ + 0(Ar- V2) (6 6) 

^U7r|e|i/2+ 40F (2Ar)i/3 + ^^'^ ^- ^^•^'^ 

We now compare (j6.5j) and (|6.6j) to the — > — oo expansion of (|6.2|) and (|6.4|) respectively. In 
relation to the latter, a straightforward calculation using the leading two terms of ()2.19|) shows 
that for 1^1 — > oo, 

roo 1/12 41 12 \ 

Using this, together with H2.19() itself, we find that ()6.5() agrees with the ^ — > — oo expansion of 
(|6.2|) for the first two terms of the 0(1) part in H6.5() . but the rational factor of ^ should be j 
for the third term. At 0{N-^^^) agreement is obtained with the first term. Similarly, we find 
that the the first term at each order in 1)6. 6|) agrees with the ^ — > — oo expansion of ()6.6() . These 
results are all consistent with the matching hypothesis. 
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The LOE and LSE 

As in the Gaussian cases, we begin by expressing the densities for the LOE and LSE in terms 
of the corresponding densities for the LUE. 
According to (|On)) and (P?7)) 

lp(X;OEA,(:E(°-i)/2e-2^^)) ~ UE^(x"e-4^^)) 
/cos2^iV,c.(^) 1 \ 1 



^V4^X(1-X) 27r{X{l- X)y/2) N' ^^'^^ 
according to (|5^ and (|2.in|) 



^ + 7-1^; OE;v(x("-i)/2e-2^-)) ~ -^p{l + — UE^(x"e-4^-)) 



(2iV)2/3^^ (2iV)2/3' >> (2Af)2/3'^^ (2A^)2/3^ 



according to (|OT]) and l(T7|) 



5o5jvW , COs2^2JV,a(^) , 1 



2(7riV)i/2x3/4(l-X)i/4 87riVX(l-X) 47riV ^^(1 - X) ' 
according to (|^ and (fmH) 



(4Af)2/3^^ (4A^)2/3' (4A^)2/3^^ (4A^)2/3^ 

-^Ai(e) ^°'Ai(.)d. + ^g±lj,Ai'(e) j^°°Ai(t)^t-^^^(Ai'(e))^ (e.io) 

Substituting H2.11() in 1)6 .71) and 1)6.9^ and expanding gives 

COS (4|C|3/2/3) (l + asin(4|C|3/2/3)) 1 ^^^^^ 

47r|e| 2vrv^ (2iV)V3 ^ • ^ 



and 



(4Ar)2/3^^ (4Ar)2/3' (4A^)2/3^^ (4A^)2/3 



^ sin(2|e|3/2/3 - 37r/4) ^ cos(4|C|3/2/3) 



2Vi|e|i/4 47r|e| 

+ 2(4iV)V3vrv1^ + + - - ^ -(4|^I'/V3)^.12) 

On the other hand let us expand (|6.8|) and (|6.in|) for ^ — > — oo. Doing so we find agreement 
with the first term at each order in H6.11|) and ()6.12p respectively, as consistent with the matching 
hypothesis. 
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6.2 Microscopic delta functions 

The results of Sections 4 and 5 tell us the asymptotic expansion of the global density, and the 
soft edge density. Here we would like to relate these expansions to the result (|1.4p and its 
Laguerre analogue. 

Consider first the Gaussian cases. For |^| large but otherwise arbitrary, we write 



p(x;ME7v(e"^'^'/2))a(x)dx 

) 

[ + [ )p(X;ME^(e-^^"'))a(X)(iX 
Jri Jr2 ' ^ 



iRi 
1 / 

+ 2iV^( J '°(l + y/2^^'^';MEjv(e-^^"'))a(l + y/2iV2/3)dy 



p(l - y/2iV2/3. ME^(e-^^^'))a(-l - y/2N^^^) dy) (6.13) 

-oo 

where MEjv = OEat, UE^v, SEat respectively, Ri = (-1 + |^|/2iV2/3^ ^ _ |^|/2Ar2/3) ^nd R2 = 
(—00, oo)\i2i, d{x) = a{Nx) and we have used the fact that p is even. Because 

p-f*(y;ME;v(e-^^"')) := Y.^^,^p{l + y /2N^/^;MEN{e~^''-')) (6.14) 

is an 0(1) quantity, we see that to leading order the second and third integrals in (|6.13|) con- 
tribute 

/oo 
p'°'\y-ymN{e-^''^"))dy. (6.15) 

However, in relation to the first integral on the right hand side of (|6.13|) . we know that terms 
which are different order in N in the asymptotic expansion of p(X;ME7v(e-^^^ )) can con- 
tribute at the same order upon the substitution 1)2.11(1 . Unlike the situation at the edge, the 
asymptotic expansion of this integral does not therefore correspond directly to the asymptotic 
of the integrand, leaving us without a method of analysis. Nonetheless some insight into the mi- 
croscopic origin of the delta functions in (|1.4|1 can be obtained as a consequence of the functional 
form of (|Fm|) for /? = 1,4. 

For P = 2 we read off from ()2.9() 

/°^*(2/;UE^(e-2A^-')) ^ (Ai'(y))2 _ y(Ai(y))2, 



while (|2.18() tells us that the leading y — > —00 behavior is 2-^/|y|/7r so ()6.15|) diverges for |^| — > 00. 
Because of the result p.4j) it must be that this is exactly canceled by a contribution from the 
bulk, and thus the edge terms (|6.15() cancel. 

For P = 1 and 4 we observe from (|5.1() and (|5.4() that /9''°^*(y; UEjv(e~^^^^) appears as 
an additive component in the scaled soft edge density, together with a further term in both 
cases. The further term has the property that it is integrable for y ^ —00. Thinking then of 
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the decomposition ()6.13|1 for |^| large, and ignoring the contribution from the non-integrable 
additive component, this then suggests that 

r p{x; 0Eiv(e-^''/2))a(x) dx ^ N C pY^{X)d{X) dX - — C "^^'^ dX 
J-oo J-1 27rJ_i^l_x2 

+(a(l)+a(-l))-y Ai(y)(l-y Ai(t)dt)(iy (6.16) 

y"°° p{x; SEjv(e~2^''))a(a;) ~ iV p^Jv{X)a{X) dX + ^ ^ "^^^^ dX 

-(a(l) + 5(-l))i r Ai(y)( r Aiit)dt)dy. (6.17) 

Here the bulk contributions are the leading two non-oscillatory terms exhibited in (|4.Hj) and (|4.17j) 
respectively, and the edge contributions are the leading terms in (|5.1|) and (|5.4j) respectively, 
with the component corresponding to the (3 = 2 edge deleted. Since 

A,fe)(l-y^ AiWd.) = --(l-y^ Am<it) (6.18) 

/■oo I d / f'°° \2 

Ai(y) y Ai(t)(it = -2^(y Ai(t)a!tj (6.19) 

we see from Ai(t) dt = 1 that the final integrals in l|6.16jl and (|6.17|) are both equal to 1/2, 
thus reclaiming (|1.4)) . 

We consider now the Laguerre analogue of H1.4|l . Let us introduce the so called chiral matrix 
ensembles chMF,N{g{x)) by the eigenvalue p.d.f. 

1 ^ 

1=1 l<j<k<N 

The simple change of variables shows 

^p(X; chME^(x2('^+(2-/3)/4)e-27V/Jx2^) ^ ^^(X^; ME^(x"e-2^'3")), (6.21) 

and thus that the Laguerre ensembles can be viewed as a Gaussian version of the chiral ensemble, 
generalized by the factor x2("+(2-/3)/4). We see from (jOH, (glSni), (E^j) and (|OTll that 

p(X;chOE^(x'^/2e-2^^^)) ~ 2/,w(X) + 
p(X;chUE^(x'^e-^^"')) ~ 2pw(^) + 



vriVVl -^2 

p(X;chSE^(x2«e-«^"^^)) ~ 2pw(^) + 



2a-8Nx\-, _ n„ ^ x-A , 1/4 



ttNVI - X^ 

where only non-oscillatory terms are included. In the case a = these expansions are precisely 
the same as for the corresponding Gaussian ensembles (the leading term in the chiral ensembles 
is 2py/{X) rather than p^{X) because X G (0, 1) rather than (—1, 1)). 
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At the soft edge X = 1 of the chiral ensembles, the fact that the scaled densities are the 
same as for the Gaussian ensembles tells us that there is a contribution 

to the smoothed density. And, although not the focus of attention of the present work, for the 
Laguerre and thus chiral ensembles there is an edge effect at X = (the so called hard edge [S]) 
which one expects to give a microscopic contribution 

2N ^ ' 

(see the next subsection) where f{X)6^ {X) dX = /(O). Consequently we expect the La- 
guerre analogue of ()1.4() to be 

p(X;chME,v(:E'^"/'e-2/5^"')) 

~ 2pw(^) + " - —6+(X) + —(l--)(6(X -I)- - , ^ ) .(6.22) 

6.3 Macroscopic balance 

In this final subsection, we will show that the results H1.4|l . 1)6. 22[) are consistent with macroscopic 
considerations. 

In a one-component log-gas, to leading order in the electrostatic potential created by the 
particle density cr(x) must balance the background potential u{x), and thus the equation 

u{x) + C = J a (y) log \x — y\dy, x€(— c, c) (6.23) 

where C is some constant, must hold. This is to be solved subject to the particle conservation 
constraint ^ 

^' a{y)dy = l. (6.24) 



For example, with u{x) = x^, (|6.23|) and H6.24|l are satisfied with 

<y{y) = Pw(y) (6.25) 

(see e.g. jSj). To the next order in A^, fluctuations in the particle density create an electric field 
and thus a force density which must be balanced for the system to be stable. The balancing 
force is provided by the gradient of the pressure fluctuation, and leads to the refinement of 1)6. 23p 

m 

u{x)+C= j a{y)\og\x -y\dy + —y-- -j\oga{x), xG(-c,c) (6.26) 
which again is subject to 1)6.24(1 . With u{x) = x^, setting 

<y{y) = Pv^{y) + ^ (6.27) 
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we see that 

C = y Ky)^og\x -y\dy + (^^- log pw{x), xG(-c,c) (6.28) 
which must be solved subject to the constraint 

y ^(y) dy = 0. (6.29) 

Differentiating 1)6. 28|) gives 

where PV denotes the principal value. Making use of the fact that 

PV / — j=L=dy = 0, xG(-l,l) 

7-1 X-y ^l-y2 

(see e.g. ^^), we see that (|6.30p . (|6.29() is solved by 



/3 2A2V V vryr^^ 

in agreement with (|1.4j) . 

The chiral ensemble (|6.20p is a log-gas confined to x > with image charges in x < 0. At 
leading order in N balance of the electric field requires 

u{x) + C= I a{y)\og\x^ -y'^\dy, xG(-c,c), (6.31) 
Jo 

subject to the constraint 

/ cr(y) dy = 0. 

But (|6.31l) can be written 

u{x) + C = j a {y) \og\x — y\dy 

so we have essentially the previous situation. In particular, with u{x) = 2x^, this shows 1)6. 3ip 
is satisfied with 

a{y) = 2p^{y). (6.32) 

To next order in N the chiral ensembles in ()6.22|1 have u{x) = 2x'^ — ^ logx, and the general- 
ization of 1)6. 28(1 reads 

-|logx + C = J ^fi{y) log \x -y\dy+(^^- logpw(2;)- 
With /i(y) I— > fi{y) — ^6{x) this is in fact identical to 1)6. 28() and we thus reclaim 1)6. 22() . 
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